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Abstract 

We use vector flavor symmetry to relate form factors of isospin changing operators to isovector 
form factors. Flavor twisted boundary conditions in lattice QCD thus allow isovector form factors 
of twist-two operators, e.g, to be computed at continuous values of the momentum transfer. These 
twisted boundary conditions, moreover, are implemented only in the valence sector. Effects of the 
finite volume must be addressed to extract isovector moments and radii at zero lattice momentum. 
As an example, we use chiral perturbation theory to assess the volume effects in extracting the 
isovector magnetic moment of the nucleon from simulations with twisted boundary conditions. 
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I. INTRODUCTION 



Numerically simulating quantum chromo dynamics (QCD) on a spacetime lattice enables 
the study of hadrons from first principles. Hadronic properties computed from lattice QCD, 
however, suflFer from a number of artifacts due to the approximations involved in solving the 
theory numerically. These artifacts include volume effects, lattice spacing errors, the use of 
unphysically large quark masses, and the partially quenched approximation. There has been 
considerable recent eflFort to address the reduction of these systematic errors using eflFective 
field theory (EFT) methods, see, e.g. [1]. With improved numerical algorithms and enlarged 
computing resources, we are entering a period where lattice data in conjunction with EFTs 
will provide physical predictions for QCD. 

Another restriction in lattice simulations is the available momentum. With periodic 
boundary conditions, the momentum and hence momentum transfers are quantized. On 
current dynamical lattices the lowest non-zero momentum mode is ^ 500 MeV. Chiral EFTs 
predict the momentum transfer dependence of form factors. With current lattice sizes, 
however, the applicability of these theories to describe the momentum dependence is ques- 
tionable. The restriction to quantized momentum transfer leads to a further impediment: 
the extraction of radii and moments that are only accessible in the near-forward limit is 
severely limited without ad hoc models for the momentum dependence. For a lattice of 
uniform spatial size L, periodic boundary conditions yield momentum modes q = 27rn/L, 
for n G Z^. To reduce the momentum granularity one must increase the lattice volume and 
thus generate new gauge configurations — an extremely costly solution.^ 

The use of periodic boundary conditions is one of simplicity not necessity. The quark 
fields need to periodic only up to some transformation which is a symmetry of the action. 
Thus if [/ is a symmetry of the action and WU = 1^ we can choose a twisted boundary 
condition for the generic field of the form 

(l){xi + L) = U(/){xi), 

while maintaining the single valuedness of the action. Twisted boundary conditions are by 
no means new H, 0,13, and there has been renewed interest in their 

utility at producing continuous hadron momentum 0, Q, 0, 0, 0, 0, 0, [20] . 

Obtaining hadronic states at continuous values of momentum does not solve the problem 
of coarse grained sampling of form factors, as the momentum transfer will still be quantized. 
In it was shown that matrix elements of fiavor changing operators, however, can be 
accessed at continuous values of the momentum transfer of the form q = 27rn/L + 86 / 
where 50 is a continuous parameter and is equal to the diflFerence of twist angles of the 
fiavors changed. As such fiavor- changing operators have no self-contractions, the difference 
of twists is implemented only in the valence quark sector to produce momentum transfer. 
These so-called partially twisted boundary conditions eliminate the need to regenerate gauge 
configurations for each value of the twisting parameters, and allow moments and radii of 
fiavor changing operators to be extracted at zero lattice momentum (up to volume correc- 
tions ll5|). This procedure has been studied in numerical simulations: for extracting /^^ jl^ . 
and for determining K ^ n matrix elements |19|]. 



One might reason that continuum moment equations could be employed to deduce moments and radii, 
e.g., using x x J{x) to determine magnetic dipole moments, or x^J4,{x) for charge radii. On the lattice, 
however, these operators do not circumvent the restriction to quantized momentum 0. 
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In this work, we point out that the vector flavor symmetry SU (2)y of QCD relates matrix 
elements of isospin changing operators to matrix elements of isovector operators. Thus the 
form factors of the latter can be accessed at continuous momentum transfer by calculating 
flavor changing matrix elements on the lattice. This result holds for quark bilinear operators 
of arbitrary spin and Lorentz structure. Our discussion is organized as follows. In Sect. ITTl 
we derive the relations between matrix elements related by an isospin rotation. Next in 
Sect, mil we detail how these relations can be utilized on the lattice with partially twisted 
boundary conditions. The dynamical effects of the boundary conditions are discussed in 
Sect. II VI where the volume effects for the nucleon isovector magnetic moment are presented. 
We end with a brief summary in Sect. 



II. VECTOR CURRENT CONSERVATION 

The Lagrangian of two flavor QCD,^ the quark part of which is 

2 

with Q = [u^dY and m = diag(m^,m^), has an exact SU{2)y symmetry in the isospin 
limit, rUu = rrid^ that cannot be spontaneously broken. We shall work exclusively in the 
isospin limit. Let us denote the generators of SU{2)v by T^. Noether's theorem yields the 
current 

j;{x) = Q{x)T'^j,Qix), (2) 

with conserved charges 

= / dxr^x), (3) 



that are the generators of isospin rotations. 

Now consider a quark bilinear operator of the form 

0^{x) = Q{x)T^TQ{x). (4) 

The r represents any Dirac matrix and any Lorentz tensor. For example, the twist-two 
operator 

= Qi^) 7,7i,D,, . . . D,„}Q{x), (5) 

is such an 0^{x)^ etc. We leave F unspecifled because only the flavor structure is relevant 
for our discussion. It is straightforward to show that 

=z5"^^(9^(x). (6) 

Deflning the usual isospin raising and lowering operators = ± zT^, we have 

0^ix) = T[Q^,0'{x)]. (7) 



^ In this Sect., we limit our discussion to the two flavor theory for simplicity. We will explain below how 
the result generalizes to the relevant subgroup of the partially quenched theory and analogous three-flavor 
theories. 
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These relations enable us to relate isospin changing matrix elements to isovector ones. 
As an example, consider the neutron to proton matrix element of . We have 

{p\0^\n) = {p\0^\p) - {n\0^\n), (8) 

by virtue of Eq. ([7j). This result can be rewritten in various forms using the isoscalar 
combination 

0'{x) = Q{x)irQ{x), (9) 

and the fact that 

{p\0^\p) - {n\0^\n) = 0. (10) 

For example, {p \ n) = 2{p\0^\p) = {p\ur u\p) — {n\ur u\n). There is a particularly 
useful way to rewrite the above relation in Eq. (jSj) for the case of the electromagnetic current 
operator 

j;^{x) = Q{x)Q-f,Q{x), (11) 
with Q = diag(2/3, -1/3). This form is 

{p I u^, d\ n) = {p\ j;^\p) - (n I J7I n). (12) 

Of course one is not limited to baryonic matrix elements. Below we shall largely highlight 
the calculation of nucleon electromagnetic form factors as an example of how to utilize 
twisted boundary conditions. Our results, however, apply to the arbitrary quark bilinear 
operators in Eq. (^, using the relation in Eq. ([Zj) between the desired states. 



III. IMPLEMENTATION ON THE LATTICE 

The above operator relations can be utilized to study isovector form factors at continuous 
values of the momentum transfer. In this Sect., we describe how to utilize partially twisted 
boundary conditions for this purpose. Here we focus on the kinematical effects; while in 
Sect. IIVI we take up the dynamical effects at finite volume. 

To see that the relations implied by Eq. ([Zj) [or for a specific example, the relation in 
Eq. (jl2j] ] allow lattice matrix elements to be accessed at continuous momentum transfer, 
first observe that there are no operator self contractions, thus, the momentum injection 
occurs in the valence sector alone. To separate the valence and sea sectors, we use the 
partially quenched Lagrangian 

6 

c= (13) 

The six quark fields transform in the fundamental representation of the graded SU{4:\2) 
group and appear in the vector Q = (2/, rf, j, /, u, c))^ that, in addition to the u and d 
quarks, has ghost quarks u and 5, which cancel the closed valence loops, and two sea 
quarks j and /. In the isospin limit, the quark mass matrix of SU{4:\2) reads tuq = 
diag(m^, m^, rrij.mj, m^, m^). 

In a finite box, the quark fields must satisfy boundary conditions that preserve the single 
valuedness of the action. Such a choice is afforded by twisted boundary conditions of the 
form 

Q{x + Lsr) = exp (ir • e^T^) Q{x), (14) 
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FIG. 1: Contraction encountered in the neutron-to-proton correlation function Eq. ((19j) . The 
source and sink are denoted by fihed circles, the operator insertion by an open circle. Single lines 
represent the propagators of twisted d-quark fields, while double lines represent the propagators of 
twisted u-quaik fields. 

where is a unit vector in the r*^ spatial direction and the block diagonal form of the 
supermatrices is 

T^ = diag(r<5,o,r5). (15) 

Here are the elements of the Cartan subalgebra of U{2).^ As a consequence of Eq. f[THj) . 
the sea quarks are periodic, and hence the twist angles 0^ can be altered without generating 
new gauge configurations. 

Defining new quark fields as Q{x) = {x)Q{x)^ where V{x) = exp{i0^ • x T^/ L)^ we can 
write the partially quenched Lagrangian as 

C= 5]Q^ (^+mg) 'g,, (16) 

where all Q fields satisfy periodic boundary conditions, and the effect of twisting has the 
form of a U{1) gauge field: = + iB^, with = {Q.O"^ T^/ L). For convenience 
we treat the twisting in the flavor basis: O^T^ = diag(0^, 0^, 0, 0, 0^, 0^), and similarly 
for Bj^ = diag(5^, 5^, 0, 0, 5^, 5^). The constant field Bj^ acts as flavor-dependent field 

momentum. Meson and baryon fields formed from these Q quark fields also acquire flavor- 
dependent momentum via the background U{1) field namely, 

D^t = d^t + i[B^,t], (17) 

for the mesons ^^1, and 

[V,B{x)f^' = d,&^\x) + i{B\^ + Bf^ + Bl)&^\x), (18) 



for the baryons [Jj 

In analogy with QCD, we assume that the infinite volume theory described by Eq. fITHj) 
will undergo spontaneous symmetry breaking of the form SU{A\2)l®SU{4\2)r SU{4\2)y. 



^ Any generator of the U(2) algebra can actually be chosen for the twists because we work in the isospin 
limit. In this case, the boundary can change i^-quarks into d-quarks and vice versa. Electric charge 
conservation requires that the electromagnetic current couple to the background field B^. Additionally 
the valence quark propagators become flavor non-diagonal, but these flavor rotations could be used to 
inject momentum transfer. For ease of applicability in lattice simulations, however, we have chosen to 
work in the Cartan subalgebra, where flavors do not rotate at the boundary, but isospin symmetry implies 
operator relations between flavor diagonal and flavor rotated currents. 



The remaining vector symmetry is explicitly broken by the quark mass difference rrij ^ rriu 
from SU{A\2)v 5[/(2|2)y ® SU{2)v The valence generators of the graded SU{2\2)v 
symmetry again lead to the operator relations in Eq. ([7j). Furthermore matrix element 
relations, e.g. that in Eq. f[T^ . remain valid because the external states contain only valence 
quarks. 

Having related isovector matrix elements to isospin changing matrix elements in partially 
twisted QCD, we briefly recall how the latter can be determined at continuous momentum 
transfer from lattice correlators. To implement the twisted boundary conditions one uses 
interpolating flelds for the hadrons built of the Q flelds that are periodic but coupled to 
the background fleld S^. For example, let us denote neutron and proton interpolating flelds 
constructed in this way as M{x^ t) and V{x^ t). One then calculates the correlation function 

= 5^e-^^•"'(0|^(^I^,^)a+(^I^^O^(0,0)|0), (19) 

cc,cc' 

where is the valence isospin raising operator 

d^{x,t) = u{x,t)td{x,t), (20) 

and the Dirac and Lorentz structure f is arbitrary. The tilde represents that any derivatives 
Djj^ in r appear as Djj^ in F. A typical contraction contributing to this correlation function 
is depicted in Fig. The kinematic effect due to twisting is uncovered by expressing the 
correlation function calculated on the lattice in terms of operators Af and V built from the 
twisted quark flelds Q. We have 

CM = {vio,t)d,it')JfiP,o)) 

= {r{Br,t)O^{t'W{P + B^,0)), (21) 

where P = 27rn/L is the lattice momentum of the neutron, Bj^ = (0^ + 20^)/L, and 
Bjy = (20^ + 0^)1 L are the fleld momenta of the neutron and proton, respectively. Notice 
the momentum transfer q = {6^ — 0^ — 2i\ri)jL can be varied continuously. 

Some flnal points are in order. Firstly, isoscalar quantities (which are notoriously difflcult 
to calculate on the lattice) of course cannot be deduced from these techniques. Next, the 
results trivially extend to the S'[/(3) and SIJ{^\?>) flavor groups because they contain the 
valence SU{2^y isospin subgroup. There are additional matrix element relations in the limit 
of an exact valence S'[/(3)y symmetry; but because this symmetry is badly broken by the 
light-quark mass difference, we shall not investigate these relations. Finally, at flnite volume, 
the presence of the By, fleld in the Lagrangian breaks the SU{2\2)y symmetry and leads to 
modiflcation of our results. These modiflcations can be addressed systematically in chiral 
perturbation theory (xPT). 

IV. NUCLEON ISOVECTOR MAGNETIC MOMENT 

In this Sect., we consider speciflcally the isovector magnetic form factor of the nucleon. 
We flrst show the limitations in extracting the magnetic moment from lattice data at the 
smallest available lattice momentum transfer. To show this is remedied by twisted boundary 
conditions, we calculate the isovector magnetic form factor in heavy-baryon xPT at flnite 
volume with partially twisted boundary conditions. 
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A. Momentum Extrapolation 

In principle the momentum dependence of the nucleon form factors predicted from heavy 
baryon xPT 0, [i^, [H, 2^ can be used to extrapolate lattice data down to zero momentum 



transfer in a model-independent way. In xPT, the isovector magnetic form factor has the 
one-loop behavior [H, Q [H, ^ 

F,{q') = 2f,^ - ^ [\xm^P^{x^q')-^^^ [\x F[m^P^{x,q'), A], (22) 
^^J Jo / Jo 

where P^rix^ q^) = ^1 + x(l — x)g2/m^, and the non- analytic function F(m, S) is given by 



F{m, 5) = -5 log g + VS^^^ log ^^-^^ ■ (23) 

Using estimates of the low-energy constants: gA = 1-25, Is'aatI = 1-5, M = 0.94 GeV, 
A = 0.29 GeV, and /ii = 3.38; we plot the deviation from linearity of the one-loop result 
for the F2 form factor in Fig. [21 for various values of the pion mass. This is done using the 
function AF2(g^) defined by 

The plot shows considerable deviation from linearity at the physical pion mass. For lattice 
pion masses around 0.25 to 0.35 GeV, there is a ~ 5 — 10% deviation from linear behavior. 
As the lattice pion masses are brought down at fixed lattice spacing (so that the minimal 

remains fixed), there is a clear trend toward non-linear behavior in g^. 

We must keep in mind that the results plotted in the Figure receive ^ 30% corrections 
for pion masses around 0.35 GeV due to higher-order terms in the chiral expansion that scale 
as m^/A^ [27, 28]. Such corrections are independent of g^. Beyond ^ 0.1 GeV^ there are 
^ 30% corrections from recoil q^ /M^ terms in loop graphs that contribute to the deviation 
from linearity [2^ . The Figure shows that non-linear q^ behavior is to be expected from the 
magnetic form factor for moderate pion masses. The EFT, however, is at a loss to describe 
the lattice data at the minimal lattice momentum transfer available on current dynamical 
lattices: q^ = 0.25 GeV ^ lies far off the axis in the plot, certainly too far to utilize the EFT. 
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B. Twisted Boundary Conditions and Finite Volume Modifications 

From our discussion in Sect. lllll we know that the restriction to discrete lattice momentum 
can be circumvented by using partiaUy twisted boundary conditions and calculating matrix 
elements of the isospin raising operator, see Eq. f[T^ . We demonstrate that this is the 
case by calculating the isovector magnetic form factor using partially twisted baryon xPT. 
Additionally by using this theory at finite volume, we can deduce the dynamical effects due 
to the boundary conditions. These effects are systematic errors that must be removed to 
determine the form factor. 



1. Partially Twisted and Partially Quenched xPT 

In the meson sector of partially quenched xPT (PQxPT) jsO, 31, [H, [s^, 01, the coset 
field S, which satisfies twisted boundary conditions, can be traded in for the field S defined 
by = V\x)Ti{x)V {x) ^ which is periodic at the boundary In terms of this field, 

the Lagrangian of PQxPT appears as 

C = ystr (p^tb^t^^ - A str (^m^S + S^mg) . (25) 

The action of the covariant derivative is specified in Eq. f[T7j) . 

To include baryons into PQxPT, one uses rank three fiavor tensors jl^l- The spin-| 
baryons are described by the 70-dimensional supermultiplet B^^^ ^ while the spin-| baryons 
are described by the 44-dimensional supermultiplet T^^^ [36] . The baryon fiavor tensors we 

use, however, are twisted at the boundary of the lattice. Thus we define new tensors 
and T^^^ both having the form ^| 

4fc(a;) = V\{x)Vy{x)Vl{x)B,,k{x). (26) 

These baryon fields satisfy periodic boundary conditions and their free and leading-order 
interaction Lagrangian has been given in fl7l |. 

In partially quenched QCD, the isovector vector current is defined by J^{x) = 
Q{x) T^^ij^Q{x). The choice of supermatrices is not unique js?'], even when one imposes 
the condition str = 0. One should choose a form of the supermatrices that maintains the 
cancellation of valence and ghost quark loops with an operator insertion jsl, . For the 
fiavor changing contributions we consider below, however, these operator self-contractions 
automatically vanish. For our calculation we require the action of in only the valence sec- 
tor, and specify the upper 2x2 block of to be the usual isospin generators T^. Henceforth 
we restrict our attention to the operator = + iJfi^- 

In calculating the isovector magnetic moment of the nucleon at next-to-leading order in 
the chiral expansion, there are local operators that contribute at tree level. In xPT, the 
leading isovector current operator is 



= j^d^ {Na,^T^N) . (27) 



M 

In partially twisted PQxPT, there are two terms in the leading isovector current operator 



SJt = ^D,. U (Ba^,BT+) + //^ (Bu^.T+b)] . (28) 
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FIG. 3: One-loop contributions to the isospin transition matrix elements in partially twisted heavy 
baryon xPT. Nucleons (deltas) are represented by single (double) lines, while mesons are repre- 
sented by dashed lines. The cross is the hairpin interaction, and the wiggly line shows the insertion 
of the isospin raising operator. 

The contribution of these operators at tree level is proportional to the linear combination 
I Ma — which is identical to the xPT low-energy constant /ii as can be demonstrated 
by matching [s^. 

2. Partially Twisted Isovector Magnetic Form Factor 

In the infinite volume limit, the isovector magnetic moment can be extracted from the 
matrix element 

{p{q) i \4\n{0) T) = ^F,{q'), (29) 

in the case where q = (0,g, 0). On the lattice, we can take both the source and sink to 
be at zero momentum, so that q = 0. Momentum transfer can then be induced by giving 
the valence up and down quarks different twist angles. For simplicity we choose = 
and = (0, 5, 0). Calculation of the infinite volume isovector form factor then proceeds 
similarly to that above in Sect. IIVXI In partially twisted baryon xPT, there are additional 
diagrams that contribute involving the hairpin interaction. These diagrams are shown in 
Fig. El The sum of all hairpin diagrams, however, vanishes. With ruj = rriu^^ the infinite 
volume contributions from the diagrams in the Figure are identical to F2{q^) in Eq. ((22|) 
under the simple replacement q —> B. This is the kinematic effect we expect from raising 
isospin with a twisted 2i-quark. 

Dynamical effects due to twisted boundary conditions arise from the propagation of the 
light Goldstone modes to the boundary. The sensitivity of these modes to the boundary 
conditions must be taken into account and can be done so in a model-independent way using 
baryon xPT in finite volume. The modification to the effective theory is straightforward. 



^ When rrij ^ m^, the results are identical to the partially quenched version of F2{q^)^ the form of which 
can be inferred from expressions in ^36i] . 



9 



FIG. 4: Additional one-loop contributions to the isospin transition matrix elements. Diagram 
elements are the same as in Fig. IHl These diagrams vanish for periodic fields in a finite volume. 



In effect, we replace the integrals over loop momenta with sums over the allowed^ modes 
k = 27Tn/L. The twisting is already taken into account in the effective theory by the U{1) 
gauge covariant derivative Z)^ above. The Poisson re-summation formula then allows us to 
cast these sums into the infinite volume result plus the finite volume modification. 

There are, however, further contributions from using the partially twisted chiral theory 
in a box.^ Diagrams that ordinarily vanish in infinite volume can now make contributions 
in a finite volume. This is the case for the diagrams depicted in Fig. 01 Quite interestingly 
these diagrams are only non-vanishing in a finite volume with twisted boundary conditions. 
When the twisting parameters vanish, so too does this finite volume effect. This can easily 
be explained. In infinite volume the diagrams in Fig. ^ vanish due to SO (4:) rotational 
invariance, while in a periodic finite volume they vanish due to invariance under lattice 
rotations. Lattice rotational invariance is broken in the direction of the twisted boundary 
conditions, hence the diagrams make non- vanishing contributions. Notice that the hairpin 
diagrams each vanish because the flavor-neutral mesons are additionally neutral under S^. 

Combining the infinite volume and finite volume results with a twisted valence 2i-quark 
and specifying the case mj = rriu^ we arrive at 

{p{0) i \4\n{0) T) = ^ U{B^) - ^f^/C.(m., By^ 0) - ||?^/C.(m., By^ A) 



3M 



47rV 



^ dx i^g\C33[m^Pn{x, B^),xBy,0] 



+ lglNJC^s3[mMx.B^),xBy,A] \ ), (30) 



where F2{B'^) is given in Eq. (|22I)- The effects of the finite volume are encoded in the 
functions lC2{m^ A) and Cssi'^i B^/S)^ which are defined in the Appendix. In the limit 
that S ^ 0, the result above accordingly vanishes: B functions as a momentum transfer 



^ As is customary we treat the length of the time direction T ^ so that for = (/cq, /c), we can take ko 

to be continuous and. A; quantized as above. 
^ As with the pions the proton and neutron are no longer degenerate due to finite volume effects. The 

volume induced isospin splittings are largest for small pion masses and grow with 0. In a 2.5 fm box at 

6> = TT, and at the physical pion mass, the splittings are ~ 15% for the pions and ~ 5% for the nucleons. 

When the pion mass is twice as big, the splittings are ~ 1% for both and are hence neglected. 



10 




L[fm] 

FIG. 5: Comparison of finite volume effects for the isovector magnetic moment. Plotted versus L 
is the relative difference foi" a few values of the pion mass. The total finite volume effect is 

denoted by K + L, while the contribution to from d]C2/dB alone is denoted by K. 



which is necessary for the magnetic form factor to be visible. Taking the derivative 

Um^(p(0)i|J3+|n(0) t), 

we can extract the magnetic moment ^2(0) up to additive volume corrections. The correc- 
tions at S = involving the £33(7717,-, 0, A) function are identical to the finite volume results 
in [40].'' Those results, however, were derived under the assumption that a momentum ex- 
trapolation to zero had been performed, or alternately that one had employed a background 
magnetic field. Devoid of these assumptions, the current insertion method produces a second 
finite volume correction to the magnetic moment involving 



d)C2{m^,By,A) 



OB 



B=0 



To investigate the effect of the finite volume on the extraction of the isovector magnetic 
moment using twisted boundary conditions, we define the relative difference 

(31) 

In the limit 5^0, the difference is just 5l[^2(0)] = Sl[ij]^ the relative difference in the 
magnetic moment. In Fig. El we plot as a function of L for various values of the 

pion mass to contrast our results with those of Ref. [40]. We plot the total contribution 
to Sl[ij] which arises from both the dlC2/dB^ and £33 functions in Eq. (|30|) . as well as just 
the contribution from dJC2/dB. The latter is the dominant finite volume effect. Of course, 
to extract the magnetic moment, we require 9^0 and thus we investigate 6l[F2{B'^)] in 
Eq. dn^. In Fig. [HI we fix the lattice volume at 2.5 fm and plot the relative difference 



^ For comparison, we have £33(771, B = 0, A) = |3^(A) of Ref. ^40j . 
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FIG. 6: Finite volume effects for extracting the isovector magnetic moment. Plotted versus 9 is 
the relative difference 5l[F2{B'^)\ for a few values of the pion mass. The momentum transfer from 
twisting is B = 6 X 0.079 GeV 

5l[F2{B'^)] as a function of 9 for a few values of the pion mass. We see that the effect of the 
finite volume decreases with 9. Generally the volume effects with momentum transfer (here 
the momentum transfer q = B = 9/L) are smaller than those at zero momentum transfer.^ 
In [39] this was anticipated due to the loop pion mass appearing as + x{l — x)q^ > m^. 
There is also a simple physical interpretation for this effect: with a space-like momentum 
transfer the correlation function is being probed on distances of order ^1/ ^/q^. As q'^ 
increases, the resolving power of the virtual probe diminishes the volume effect. 

While twisted boundary conditions have introduced systematic error into the determi- 
nation of the isovector moment, we stress that this is a controlled error. In a fixed box 
size, we need only know the low-energy constants qa^ Qan^ A, and to remove this effect. 
Without twisting, one must rely on phenomenological or other fitting functions which in- 
troduce uncontrolled error. With twisting, however, the lattice practitioner can approach 
the calculation of the isovector magnetic moment from various ways in conjunction with the 
EFT. For example, one can choose values of ^ 0.1 to eliminate the need for a momentum 
extrapolation, cf Fig. [21 One must then use the EFT to remove the ^ 25% volume effects. 
Alternately one could choose 9 < tt to minimize the effect of the finite volume. For these 
values of 9 one then uses the momentum dependence predicted by the EFT, cf Eq. (|22|), 
to obtain the magnetic moment. Another virtue to the EFT approach is that higher-order 
corrections to the quark mass, volume, and momentum-transfer dependence can be calcu- 
lated to aid in the extrapolation. Although we have focused on the magnetic moment, the 
expressions we have derived here are also relevant for the isovector magnetic radius. 

As a final point, an assumption inherent in our discussion is that no multi-particle 
thresholds are reached. For large enough momenta, multi-particle thresholds are inevitably 
reached. Volume corrections are no longer exponentially suppressed in asymptotic volumes; 
they become power law and will likely dwarf any signal. For the extraction of moments and 



The actual behavior with respect to momentum transfer is damped oscillatory. The oscillations arise from 
the kinetic term {k + xq)'^; but, in an L = 2.5 f m box, they set in beyond the reach of the effective theory. 
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radii of stable particles near zero momentum transfer, however, we are safely away from the 
multi-particle continuum. 



V. SUMMARY 

We have related matrix elements of arbitrary quark bilinear operators that change isospin 
to isovector combinations of those that do not change isospin. Using flavor twisted boundary 
conditions on the valence quark flelds, form factors of these isovector matrix elements can 
then be deduced from the lattice at continuous values of the momentum transfer. Isovector 
moments and radii can thus be determined from simulations at zero lattice momentum. 
Using twisted boundary conditions on the quark flelds dramatically modifles the effect of 
the flnite volume, even away from multi-particle cuts. This systematic effect can be handled 
with EFTs. The determination of isovector moments and radii without any model-dependent 
assumptions about the momentum dependence is thus well within reach of current resources. 
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FINITE VOLUME FUNCTIONS 

In this Appendix, we deflne and evaluate the flnite volume functions contributing to the 
isovector magnetic form factor. These functions can be related to a more basic sum that is 
ubiquitously encountered in these calculations 

^n{B, P^) = J^Yl \(h. ^ , ^2]a ~ [ 79^171^X^^X7^' ^^^^ 



^3 ^ [(fc + By + /?2]- J (2^)3 [(fc + By + 

which can be cast into an exponentially convergent form involving elliptic theta func- 
tions In the text, the function /C2(m, A) is deflned by 



/C,(m, A) = ^ ^ rfA X: [(fe+|)tf^2p/. ^ (33) 



with = + 2AA + m^. Using the expression for T1/2 in terms of elliptic theta functions, 
the integral over A can then be expressed in terms of the Erfc(x) function j4]J. The flnal 
result appears as the one-dimensional integral 



IC2im,By,A) = c^r r-^/V^^'"'"') Erfc(AV^K3 (^^,e-^^ ^3 (o,e-^)', 

(34) 
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with i?3(g, z) as the Jacobi eUiptic theta function of the third kind. Lastly the function 
£33(777,, B,A) is given by 



poo -1 

£33(m,5y,A) = SttM dX —Y1 

J^ u 



[(fc + B)2 + /?2]5/2 



r dk 

J M 



^^^-3/2gr(A2-m2) Erfc(AV^) 



)3 [(fc + B)2+/32]5/2 



0,e-^ 



-, e 4t 



8r 



+ ^^^3 (0, e"^ ) ( 0, e"^ ) i9 



-, e 4t 



(35) 
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